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1. Introduction 



The discovery of the Higgs boson and the measurement of its interactions with massive quarks 
and vector bosons represent a central goal of the Large Hadron Collider (LHC). In the mass range 
1 14GeV < Mh i$ 130GeV, associated ttH production provides the opportunity to observe the Higgs 
boson in the H — > bb decay channel and to measure the top-quark Yukawa coupling. However, the 
extraction of the ttH(H — > bb) signal from its large QCD backgrounds, pp — > ttbb and ttjj, repre- 
sents a serious challenge. The selection strategies elaborated by ATLAS and CMS [0, |], which 
assume 30 fb -1 and 60 fb -1 , respectively, anticipate a statistical significance around 2a (ignoring 
systematic uncertainties) and a signal-to-background ratio as low as 1/10. This calls for better than 
10% precision in the background description, a very demanding requirement both from the exper- 
imental and theoretical point of view. Very recently, a novel selection strategy based on highly 
boosted Higgs bosons has opened new and very promising perspectives [fsj. This approach might 
increase the signal-to-background ratio beyond 1/3. 

The calculation of the NLO QCD corrections to the irreducible ttbb background, first pre- 
sented in [|], |5|] and subsequently confirmed in [0], constitutes another important step towards the 
observability of ttH(H — > bb) at the LHC. These NLO predictions are mandatory in order to reduce 
the huge scale uncertainty of the lowest-order (LO) ttbb cross section, which can vary up to a factor 
four if the QCD scales are identified with different kinematic parameters []7p. Previous results for 
five-particle processes that feature a signature similar to ttbb indicate that setting the renormaliza- 
tion and factorization scales equal to half the threshold energy, /Ir f = £V/2, is a reasonable scale 
choice. At this scale the NLO QCD corrections to pp ->■ ttH (K ~ 1.2) [g], pp -> ttj (K ~1.1) [§], 
and pp — >• ttZ (K ~ 1.35) JIfj[], are fairly moderate. This motivated experimental groups to adopt the 
scale jUr p = E^/2 = m t + m b ^ /2 for the LO simulation of the ttbb background [jlj]. However, at this 
scale the NLO corrections to pp — > ttbb turn out to be unexpectedly large (K ~ 1.8) []5| Hj. As we 
argue, a reliable perturbative description of ttbb production requires a different scale choice [11]. 

The calculation of the NLO corrections to pp — > ttbb constitutes also an important techni- 
cal benchmark. The description of many-particle processes at NLO plays an central role for the 
LHC physics programme, and the technical challenges raised by these calculations have triggered 
an impressive amount of conceptual and technical developments. Within the last few months, 
this progress has lead to the first NLO results for six-particle processes at the LHC, namely for 
pp — > ttbb [^ ||], the leading- [12] and the full-colour contributions [|l3|] to pp — > Wjjj, and for the 
qq contribution to pp — > bbbb [14]. 

To compute the virtual corrections to ttbb production we employ explicit diagrammatic rep- 
resentations of the one-loop amplitudes and numerical reduction of tensor integrals [15, 16]. The 
factorization of colour matrices, the algebraic reduction of helicity structures, and the systematic 
recycling of a multitude of common subexpressions — both inside individual diagrams and in tensor 
integrals of different diagrams that share common sub-topologies — strongly mitigate the factorial 
complexity that is inherent in Feynman diagrams and lead to a remarkably high CPU efficiency. 
The real corrections are handled with the dipole subtraction method Jl7|]. Our results have been 



confirmed with the HELAC-1LOOP implementation of the OPP method [JTS|, [19], [20p within the 
statistical Monte Carlo error of 0.2% [Hi. 



2 



NLO QCD corrections to ttbb production at the LHC 



Stefano Pozzorini 



2. Description of the calculation 

In NLO QCD, hadronic ttbb production involves the 2—^4 partonic channels qq — > ttbb 
(7 trees and 188 loop diagrams) and gg — > ttbb (36 trees and 1003 loop diagrams). The 2—^5 
bremsstrahlung contributions comprise the crossing-symmetric channels qq — > ttbbg, gg — > ttbbg, 
and gq — > ttbbg (64 diagrams each), and the partonic process gg — > ttbbg (341 diagrams). Each 
contribution has been worked out twice and independently, resulting in two completely indepen- 
dent computer codes. The treatment of the qq- and gluon-induced reactions are described in [ffl] 



and [11], respectively. Here we mainly focus on the virtual corrections in the gg channel. 



Feynman diagrams are generated with two independent version of FeynArts [21, 22 ] and 



handled with two in-house MATHEMATICA programs that perform algebraic manipulations and 



generate FORTRAN77 code fully automatically. One of the two programs relies on FormCalc [23] 
for preliminary algebraic manipulations. The interference of the one-loop and LO matrix elements, 
summed over colours and helicities, is computed on a diagram-by-diagram basis, 



££^ 1 - 1o °p)(^ l °)) =£ 

col hel r 



col hel 



(2.1) 



Individual loop diagrams (r) are evaluated by separate numerical routines and summed explicitly. 
The cost related to the large number of diagrams is compensated by the possibility to perform 



colour sums very efficiently thanks to colour factorization [see (|2.2|)]. Individual (sub)diagrams 
consist of a single colour-stripped amplitude multiplied by a simple colour structure, 1 which 
is easily reduced to a compact colour basis {^}. The LO amplitude is handled as a vector in 
colour space, and colour sums are encoded once and for all in a colour-interference matrix 

sen = facPvA , ^ (L0) = £^ (LO V,, hi = J>*j\ (2-2) 

\ k / 1 col 

These ingredients yield colour-summed results by means of a single evaluation of the colour- 
stripped amplitude &f( T > of each (sub)diagram. Tensor integrals with N propagators and P Lorentz 
indices are expressed in terms of totally symmetric covariant structures {g...gp ... p}jl'"fp involv- 
ing external momenta pi,... ,pn-i an d g^ v in D = 4 — 2e dimensions Jl6|], 



The gg channel involves tensor integrals up to rank P = 4. The coefficients , p are related to 
scalar integrals by means of numerical algorithms that avoid instabilities from inverse Gram deter- 
minants and other spurious singularities [|^, |l6| ]. The tensor rank and the number of propagators of 
integrals with Af > 4 are simultaneously reduced without introducing inverse Gram determinants. 
Tensor integrals with ,/V = 4, 3 are handled with the Passarino-Veltman algorithm as long as no 



'More precisely, each quartic gluon coupling generates three independent colour structures that are handled as 
separate subdiagrams. However, most diagrams do not involve quartic couplings, and their colour structure factorizes 
completely. 
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small Gram determinant appears in the reduction. Otherwise, expansions about the limit of vanish- 



ing Gram determinants and possibly other kinematical determinants are applied Q16Q. The reduction 
is strongly boosted by a cache system that recycles tensor integrals among diagrams with common 
subtopologies. The gg channel involves about 350 scalar integrals, which require 10ms CPU time 
per phase-space point. 2 The calculation of all scalar and tensor integrals with and without cache 
system takes 40 ms and 200ms, respectively. Rational terms arising from ultraviolet (UV) poles of 
tensor integrals with D-dependent coefficients are automatically extracted by means of a catalogue 
of residues R N , ■ , 



(fj 



Rational terms resulting from infrared poles must be taken into account only in wave-function 
renormalization factors, since they cancel in truncated one-loop amplitudes [Q]. 

The helicity-dependent parts of all diagrams are reduced to a common basis of so-called Stan- 
dard Matrix Elements (SMEs) of the form 

An = 2,™'" %, ipi^iPi) [v( J p 3 )7 Pl • • • YpAp*)] HP5)Y Pk+l ■ ■ ■ y Pl u{j>6)] , (2-5) 
where <2m ^ 2Pl " P/ are combinations of metric tensors and external momenta. The colour-stripped 



part of each loop diagram [see (|2.2[)] yields a linear combination of SMEs and tensor integrals, 

^ (r) = £^i r) ^ m , ^"i r) =£ £ ^ ^T^ + rational parts. (2.6) 

m P h,~,jp=0 

The use of SMEs enables very efficient helicity summations. Helicity and colour sums are encoded 
in the interference of SMEs A m and colour structures ^ with the LO amplitude, 



Mhn = ££^ m ^(V (LO) ) = £/«£X(^ {L0) ) • (2.7) 



.(LO) 

col hel / hel 



This matrix, which must be computed only once per phase-space point, links the colour/helicity- 
independent form factors of each diagram to its colour/helicity-summed contribution 



col hel 

The SME-reduction starts with process-independent D-dimensional relations such as momentum 
conservation, Dirac algebra, transversality, and gauge-fixing conditions for the gluon-polarization 
vectors. Once rational terms are extracted, we further reduce SMEs with two alternative algo- 
rithms in four dimensions. The first algorithm splits each fermion chain into two contributions, 
u iPi) = La=± ®x u {Pi)i Via insertion of chiral projectors co± = (1 ± y 5 )/2. This permits to em- 
ploy various relations of type "f-y""^ (0±®y ll = (£>±® 'Y* (Q± + y"/' 3 which connect 
Dirac matrices of different fermion chains [|], |24[]. In this way a rich variety of non-trivial identities 



"All CPU times refer to a 3 GHz Intel Xeon processor. 
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are obtained that reduce the full amplitude to 502 SMEs [|Tl|]. Besides this procedure, which de- 
pends on process-specific aspects, we implemented a simple process-independent reduction based 
on a four-dimensional identity of type y^y^y^y^y^ = gWzyteyl^'ft - gWigP>iV4>fs + pe rm., 
which eliminates spinor chains with more than three Dirac matrices without introducing 75 [jnj]. 
This leads to 970 SMEs. In spite of the factor-two difference in the number of SMEs, we found that 
the numerical codes based on the two different reductions have the same — and remarkably high — 
CPU speed: about 1 80 ms per phase-space point. This unexpected result means that the obtained 
CPU performance, at least for this process, does not depend on process-dependent optimisations. 
To handle the real corrections we employed the dipole subtraction method [ 17[ 25]. The 2 — > 5 



matrix elements were generated with Madgraph pq] and checked against analytic calculations 



and in-house code based on off-shell recursions. More details are given in [11]. 

3. Predictions for the LHC 

We present results for pp — )• ttbb +X at \/s = 14TeVforra t = 172.6GeV and = 0. Collinear 
final-state partons are recombined into jets with \J Atp 2 +Ay 2 > 0.4 using a ^-algorithm. We 
impose the cuts /?T.b > 20GeV and \y\, \ < 2.5 on b jets and use the CTEQ6 set of PDFs with Np = 5 
active flavours and A^ s = 226 MeV. Further details are given in JIT)]. 

In all recent ATLAS studies of ttH(H — > bb) ||l|] the signal and its ttbb background are sim- 
ulated by setting the renormalization and factorization scales equal to half the threshold energy, 
Z?thr = 2m t + wibb- Being proportional to a s 4 , these LO predictions are extremely sensitive to the 
scale choice, and in ^ we found that at /Xr p = E^/2 the NLO corrections to pp — > ttbb are close 
to a factor of two. This enhancement is due to the fact that pp — > ttbb is a multi-scale process 
involving various scales well below E^/2. In particular, the cross section is saturated by b quarks 
with /?T,b *C m t [pT|]. In order to avoid large logarithms we advocate the use of the dynamical scale 

Mo = ra t v / Fr,bPT,b> (3.1) 

which improves the perturbative convergence and minimises NLO effects in the shape of distribu- 
tions [pi]]. In Fig. H] we show results for the kinematic region m bb - > lOOGeV, which is relevant 
for ATLAS/CMS studies of ttH(H -)• bb). The left plot displays the dependence of the LO and 
NLO cross sections with respect to scale variations /Ir = jUp = <§Ato- At the central scale we obtain 
Olo = 786.3(2)fb and Onlo = 978(3)fb. This NLO result is 2.18 times larger as compared to 
the LO cross section based on the ATLAS scale choice [O]. The scale choice (3.1) reduces the K 



factor to 1.24, and the NLO (LO) uncertainty corresponding to factor-two scale variations amounts 
to 21% (78%). The improvement with respect to [||] is evident also from the stability of the NLO 
curve in Fig. [j]a. The right plot in Fig. p] displays LO (blue) and NLO (red) scale-dependent pre- 
dictions for the m b g distribution. The results are normalized to the LO distribution at /Ir f = £io and 
the bands correspond to factor-two scale variations. The NLO predictions perfectly fit within the 
LO band and exhibit little kinematic dependence. Inspecting various other distributions we found 
similarly small NLO corrections to their shapes. In Fig. ^| we show similar plots for the kinematic 
region p Tbb > 200GeV, which permits to increase the separation between the Higgs signal and 
its ttbb background At the central scale we obtain <j LO = 451.8(2) fb and Onlo = 592(4) fb. 
The K factor (1.31), and the LO (79%) and NLO (22%) scale dependence behave similarly as for 
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Figure 2: pp — > ttbb + X cross section at the LHC for standard cuts and p T b g > 200 GeV: scale dependence 
of the LO and NLO cross section (left plot) and relative NLO corrections to the m bb distribution (right plot). 



m bb > 100 GeV. But in this case the NLO corrections are rather sensitive to m bb . In the physi- 
cally interesting region of m bb ~ lOOGeV, the shape of the m bb distribution is distorted by about 
20%. This effect tends to mimic a Higgs signal and should be carefully taken into account in the 
ttH(H -)• bb) analysis. 

4. Conclusions 

The observation of the ttH(H — > bb) signal and the direct measurement of the top-quark 
Yukawa coupling at the LHC require a very precise description of the pp — > ttbb irreducible back- 
ground. The NLO QCD corrections to ttbb production reveal that the scale choice adopted in 
previous LO studies underestimates this cross section by a factor of two. We advocate the use of a 
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dynamical scale that stabilizes the perturbative predictions reducing the K factor to about 1.2. In 
presence of standard cuts NLO effects feature small kinematic dependence. But in the regime of 
highly boosted Higgs bosons we observe significant distortions in the shape of distributions. 

The calculation is based on process-independent algebraic techniques, which reduce loop di- 
agrams to standard colour/helicity structures, and numerically stable tensor-reduction algorithms. 
The very high numerical stability and CPU efficiency of this approach are very encouraging in view 
of future NLO calculations for multi-particle processes. 
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